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We study the effects of quantum interference in impurity structures consisting of two or three 
magnetic impurities that are located on the surface of an s-wave superconductor. By using a 
self-consistent Bogoliubov-de Gennes formalism, we show that quantum interference leads to char- 
acteristic signatures not only in the local density of states (LDOS), but also in the spatial form of 
the superconducting order parameter. We demonstrate that the signatures of quantum interference 
in the LDOS are qualitatively, and to a large extent quantitatively unaffected by the suppression of 
the superconducting order parameter near impurities, which illustrates the robustness of quantum 
interference phenomena. Moreover, we show that by changing the interimpurity distance, or the im- 
purities' scattering strength, the s-wave superconductor can be tuned through a series of first order 
quantum phase transitions in which the spin polarization of its ground state changes. In contrast 
to the single impurity case, this transition is not necessarily accompanied by a vr-phase shift of the 
order parameter, and can in certain cases even lead to its enhancement. Our results demonstrate 
that the superconductor's LDOS, its spin state, and the spatial form of the superconducting order 
parameter are determined by a subtle interplay between the relative positions of the impurities and 
their scattering strength. 
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I. INTRODUCTION 

Over the last few years, quantum interference 
effects arising from scattering of conduction elec- 
trons by multiple impurities have attracted significant 
experimental^'^'^'^^Mi and theoretical^'^'-^^^-^-^'-^^^-^^'-^^ at- 
tention. For example, Manoharan et alX used a corral 
of magnetic impurities on the surface of a metallic host 
to demonstrate that quantum interference and the re- 
sulting formation of eigenmodes can be used to create 
quantum images (also called quantum mirages). Comple- 
mentary to this experiment, Derro et al? ^ using scanning 
tunneling spectroscopy (STS), reported the formation of 
novel resonance states arising from impurities in the one- 
dimensional chains of YBa2Cu306+ic. Quantum interfer- 
ence effects have also been studied in optical quantum 
corrals^, in quantum dotsi, in triangular corrals^, in fer- 
romagnetic islands^, and around molecules^. 

A number of theoretical studies have focused on quan- 
tum interference effects in metallic host systems^, as well 
as conventional^'^^'^^'^^ and unconventional^^' super- 
conductors. Superconducting (SC) host systems with 
s-wave symmetry are of particular interest for the inves- 
tigation of quantum interference effects for two reasons. 
First, a magnetic impurity induces a fermionic bound 
state in an 5- wave superconductor, which can be used as 
the object (i.e., the "quantum candle") for the formation 
of a quantum image^^, while a non-magnetic impurity 
does not. This qualitative difference between a magnetic 
and non-magnetic impurity allows one to study quan- 
tum interference effects separately from the formation 
of fermionic impurity states, in contrast to unconven- 
tional superconductorsi2ii^iiSii&. Second, quantum in- 
terference can be used to tune an 5- wave superconductor 



through a series of first order quantum phase transition 
in which the spin polarization of the superconductor's 
ground state changes^ ^'^^'^^'^ii^. Quantum interference 
phenomena are therefore not only of great fundamental 
interest, but might also possess important applications 
in the field of spin electronics^^^ and quantum informa- 
tion technology^^, as these phenomena might lead to the 
creation of new types of quantum qubits. 

In this article, we study quantum interference phe- 
nomena arising in structures consisting of two or three 
magnetic impurities that are located on the surface of 
an 5- wave superconductor. We employ a self-consistent 
Bogoliubov-de Gennes (BdG) formalism.^^ which allows 
us to study the effects of quantum interference not only 
on the formation of impurity bound states, but also on 
the spatial form of the superconducting order parame- 
ter. We obtain three important results. First, we show 
that quantum interference leads not only to character- 
istic signatures in the local density of states (LDOS) 
but also in the spatial form of the superconducting or- 
der parameter. In particular, when the distance between 
impurities is changed, the order parameter at the sites 
of the impurities (i.e., the on-site order parameter) ex- 
hibits Friedel-like oscillations. These oscillations mir- 
ror the ones in the frequency of the hybridized impurity 
bound state. Second, we find that the characteristic sig- 
natures of quantum interference in the LDOS, such as, 
e.g., the formation of bonding and antibonding bound 
states around two magnetic impurities, and the oscilla- 
tion of the bound state frequencies with interimpurity 
distance^'^^, remain qualitatively and to a large extent 
quantitatively unchanged by the suppression of the order 
parameter around impurities. This results demonstrate 
the robustness of quantum interference phenomena. The 
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physical origin of this robustness lies in the fact that, 
while the superconducting order parameter is suppressed 
by impurities, it recovers its bulk value over a length scale 
which is set by the Fermi wavelength, Ai?, and not by the 
superconducting coherence length, ^c, where for many 
s-wave superconductors typically ^ ^f- Third, we 
show that quantum interference leads to a series of first 
order quantum phase transitions, in which the spin po- 
larization of the superconductor's ground state changes. 
These transitions can be achieved by increasing the im- 
purities' scattering strength, or, which is experimentally 
more relevant, by changing the inter impurity distance. 
While this type of quantum phase transition is well un- 
derstood for a single impurity^^, we find that the phase 
transitions associated with multiple impurities exhibit 
qualitatively new features. In particular, we show that 
in contrast to the single impurity case, these transitions 
are not necessarily accompanied by a 7r-phase shift of the 
on-site order parameter, and can, in certain cases, even 
lead to its enhancement. The theoretical predictions for 
the spatial form of the superconducting order parame- 
ter discussed below are of particular importance in light 
of recent experimental progress in developing Josephson 
STS^^, which suggests that our predictions can be tested 
in the near future. 

The rest of the paper is organized as follows. In SecHll 
we briefly introduce the BdG formalism. In Sees. IIIII 
and IIVI we present our results for the case of a two and 
three impurity systems, respectively. Finally, in Sec. |V| 
we summarize our results and conclusions. 



II. BOGOLIUBOV-DE GENNES FORMALISM 

In order to treat quantum interference between elec- 
trons that are scattered by multiple impurities and a spa- 
tial variation of the superconducting order parameter on 
equal footing, we use the self-consistent Bogoliubov-de 
Gennes (BdG) formalismSi. Within this formalism, one 
solves the eigenvalue equation 




with i^jj = tij + — Here, Uj is the hopping 

integral between sites i and j, and ja is the chemical po- 
tential. Below, we consider a two-dimensional system, in 
which only the nearest-neighbor hopping integral t and 
the next-nearest neighbor hopping integral f are non- 
zero. For concreteness, we take f /t = 0.2 and /i/t = —2, 
leading to an almost circular Fermi surface (FS) with 
kp ^ 7r/2. J\ and 5'i are the magnetic scattering po- 
tential and the spin of the impurity located at site i, re- 
spectively. Below, we consider for concreteness magnetic 
impurities with spin S = 1/2. In an 5-wave superconduc- 
tor, the superconducting order parameter Ay = Ai^ij is 



local and given by 

Ai = -l^^^^i,n^i,ntanh('^|^') , (2) 

where V is the effective pairing interaction and the sum 
runs over all eigenstates of the system. Eqs.Q and Q 
are solved self-consistently in order to obtain the energy. 
En-, of all eigenstates of the system, the local supercon- 
ducting order parameter, Ai, as well as the Bogoliubov 
coefficients iii,n,'^i,n of state n at site i. Below, we take 
V/t = 2.5, such that in the clean system, the translation- 
ally invariant superconducting order parameter is given 
by Ai/t ^ 0.112 . Finally, the local density of states at 
site i is obtained via 

N{u;,i) = Y^[ulJ{L0-Er.) + vl,5{u; + En)\ . (3) 

n 

Unless otherwise stated, we consider in the following a 
system with N = x Ny sites where = 28 and 
Ny = 18. This choice of N^^Ny accounts for the spatial 
anisotropy of the impurity structures considered below. 

Note that within the BdG-formalism, one assumes that 
the impurity spin is a classical, static variable, which cor- 
responds to the limit JS = const, and S ^ oo. In a fully 
gaped s-wave superconductor, this approximation is jus- 
tified even when considering magnetic impurities with 
6* < oo, and in particular small values of S. Specifically, 
it was shown^^ that no Kondo effect occurs in a fully 
gapped s-wave superconductor for sufficiently small cou- 
pling between the magnetic impurities and the delocal- 
ized electrons. Hence, it is not necessary to consider the 
full quantum dynamics of the magnetic impurity, which 
can thus be treated as a classical, static variable. This 
finding is also supported by the experimentally measured 
LDOS around a single magnetic impurity in an 5- wave 
superconductor^^. These experiments have reported the 
existence of two peaks inside the superconducting gap, 
which are induced by the magnetic impurity. Within the 
self-consistent BdG-formalism or the non-self-consistent 
T-matrix approach^^'^^'^^ these peaks are a direct con- 
sequence of the static nature of the magnetic impurity 
and are the spectroscopic signature of a fermionic (Shiba) 
bound state. 

Finally, within the BdG-formalism, any interaction be- 
tween the magnetic impurities is only taken into account 
to the extent that it determines the angle, a, between the 
direction of the impurity spins. Within a T-matrix ap- 
proach, it was shown in Ref.^^ that for two impurities, the 
cases a = and a are qualitatively very similar^^. 
We therefore focus below on the case of impurities with 
parallel spin, i.e., a = 0. 

III. TWO MAGNETIC IMPURITIES 

Before considering the effects of two magnetic impuri- 
ties on the local electronic structure of an 5-wave super- 
conductor, we briefly review the salient features of the 
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FIG. 1: (Color online), (a) LDOS at ri,2 as a function of 
energy for two impurities located at ri,2 = (=f2,0) (interim- 
purity distance Ar = 4) with JS/t = 2 (the lattice constant, 
ao, is set to unity). Spatial intensity plot of the LDOS at (b) 
ul/t = 0.037, and (c) c^J/t = 0.079. The filled red circles 
denote the positions of the magnetic impurities. 



fermionic bound state induced by a single impurity. For 
a single magnetic impurity with S||z, particle- hole mix- 
ing in the superconducting state yields a wave-function 
of the induced bound state that possesses a particle-like 
and a hole-like component, denoted by |p, t) and j), 
respectively. Here, j, [ denote the opposite spin quantum 
numbers, Sz = ±1/2 of the two components. The spec- 
troscopic signature of the bound state are two peaks in 
the LDOS, as shown in Fig. QJa). The peak at negative 



energies arises from the particle-like component of the 
wave- function, t), while the peak at positive energies 
represents the hole-like part, \h^[). These two peaks are 
therefore also referred to as the spin-t and spin-| peaks 
in the LDOS. 

We now turn to the case of two magnetic impurities 
with parallel spins and S||z. If the impurities are in- 
finitely far apart (with inter impurity distance Ar = oo), 
the two bound states with components |p, j, i) and \ h^[^i) 
(i = 1,2) are degenerate. This degeneracy is lifted for 
Ar < oo, since an electron scattered by one of the im- 
purities is in general also scattered by the second impu- 
rity, resulting in a coupling of the induced bound states. 
In order to gain further insight into the nature of this 
coupling, we consider a toy model in which the unhy- 
bridized bound states of energy Eq are coupled by a hop- 
ping term D{Ar)^ which depends on the distance, Ar, be- 
tween the two impurities. This coupling, in turn, leads to 
the formation of even and odd bound states with energies 
£^e,o = ^0 ± D{Ar) and a particle-like component of the 
wavefunction given by |p, t)e,o = T, 1) ± |p, T, 2))/y2, 
and similarly for the hole-like component. As a result, 
the bound state peaks in the LDOS are split [by an en- 
ergy AE = Ee-Eo = 2D{Ar)], as shown in Fig.[Ha) for 
the case of two impurities, located at ri = (—2,0) and 
r2 = (2,0) (Ar = 4) with scattering strength JS/t = 2 
(in the following, we denote the energies of the spin-| 
and spin- 1 peaks of the even and odd bound state by 
cc;^'^, respectively). This splitting of the bound state 
peaks was also obtained within the non-self-consistent 
T-matrix formalismii (see also Ref.^). In Figs.^b) and 
(c) we present a spatial plot of the LDOS at and u^, 
respectively [light (dark) color indicates a large (small) 
LDOS, and the filled red circles denote the positions of 
the magnetic impurities]. Fig. ^c) clearly demonstrates 
the odd symmetry of the bound state since the LDOS 
vanishes at any point with equal distance to the two im- 
purities. 

In Fig. 13 we present the superconducting order pa- 
rameter as a function of position for the case of the two 
impurities considered in Fig. ^ As expected, we find 
that the order parameter is suppressed in the vicinity of 
the impurities, with the largest suppression at the impu- 
rity sites (indicated by arrows in Fig. [21). Note, however, 
that the order parameter recovers its bulk value within 
a few lattice spacing from the impurities. Surprisingly 
enough, even for the relatively small interimpurity dis- 
tance of Ar = 4 the superconducting order parameter 
comes close to its bulk value in the region between the 
two impurities. We therefore conclude that the supercon- 
ducting order parameter relaxes back to its bulk value on 
a length scale of a few lattice spacings. This length- 
scale Ar is much shorter than the superconducting co- 
herence length given hy = vp/A ^ 25ao, where is 
the Fermi velocity, and ao is the lattice spacing. This re- 
sult is in agreement with numerical^^ and analytical^^'^^^ 
studies of the superconducting order parameter near a 
single impurity, in which the latter identified A^ with the 
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FIG. 2: Superconducting order parameter as a function of 
position along the line connecting the two impurities at ri,2 = 
(=F2, 0) considered in Fig. Q The positions of the impurities 
are denoted by arrows. 



Fermi wavelength A^. In our case, Af ~ 4ao <C ^c- 
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FIG. 3: (Color online). LDOS for a spatially uniform super- 
conducting order parameter. The intensity plot of the LDOS 
is shown at (a) cjj/t = 0.041, and (b) c^J/t = 0.084. 



In order to ascertain the importance of the order pa- 
rameter suppression on the frequency and spatial form of 
the bound states, we computed the LDOS by constrain- 
ing the superconducting order parameter to its spatially 
homogeneous bulk value. In this case, the spin-| and 
spin- 1 peaks of the even and odd bound state are lo- 
cated at uyt = -0.042, ijl/t = 0.041, u^t = -0.086 
and oo^t = 0.084, respectively. In contrast, when a 
suppression of the superconducting order parameter is 
taken into account within the BdG-formalism, the bound 
state peaks are located at co^t = —0.038, cjj/t = 0.037, 
co^t = -0.081 and cjj/t = 0.079, respectively [see 
Fig. Q^a)]. Thus the spatial suppression of the super- 
conducting order parameter leads to a small shift of the 
bound state peaks to lower energies (a qualitatively sim- 
ilar effect was observed near a single magnetic impurity, 
see Ref.^^). In Figs. Ofa) and (b) we present a spa- 
tial intensity plot of the LDOS at ^^/t = 0.041, and 
uoyt = 0.084, respectively. A comparison of Figs. [2fa) 
and (b) with Figs. Qjb) and (c) immediately shows that 
the suppression of the superconducting order parameter 
in the vicinity of the impurities does not lead to percep- 



tible changes in the spatial LDOS pattern of the even 
and odd bound states (for better comparison, the inten- 
sity scale for all four figures is the same). Hence, we 
conclude that characteristic signatures of quantum inter- 
ference in the LDOS remain qualitatively and to a large 
extent quantitatively unchanged by the suppression of 
the order parameter around impurities. This robustness 
of quantum interference phenomena arises from the fact 
that the superconducting order parameter recovers its 
bulk value over a length scale which is typically much 
shorter than the superconducting coherence length, as 
discussed above. This conclusion also supports the valid- 
ity of the results pertaining to two impurity interference 
effects obtained in Ref.^^ within the non-self-consistent 
T-matrix approach. 

As the distance between the two impurities is changed, 
the frequency of the even and odd bound state peaks, as 
well as the splitting between them oscillates^'^^, as shown 
in Fig. 2fa). The strength of the coupling, D{Ar)^ be- 
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FIG. 4: (Color online), (a) Frequency of the even and odd 
bound state peaks, and (b) on-site superconducting order pa- 
rameter (at ri,2) as a function of inter- impurity distance Ar 
for JS/t = 2. 

tween the impurity states can be directly obtained from 
the frequency splitting between the even and odd bound 
states, AE = — Eq = 2D{Ar). The results shown in 
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Fig. 2ta) then imply that the effective couphng between 
the impurity states does not only vary in magnitude, but 
also changes sign as the distance between the impurities 
is changed. One obtains, for example, D{Ar = 2) > 0, 
while D{Ar = 4) < 0. The sign and magnitude of D{Ar) 
is determined by scattering processes involving both im- 
purities. For example, consider the wave function of an 
electron that is scattered by the impurity at ri . The sign 
and magnitude of its wave function at T2 (and thus those 
of the effective coupling) depend (a) on the distance be- 
tween the impurities due to the (A;i?r)-oscillations of the 
electronic wave-function, and (b) the scattering phase 
shift at ri , which is determined by the scattering strength 
JS. As a result of the latter, the relative splitting be- 
tween the even and odd bound states in general varies 
with JiS, as demonstrated below. Moreover, if the elec- 
tron's wave function vanishes at r2, the impurity bound 
states associated with each of the two impurities cannot 
hybridize and hence remain degenerated^. 

Complementary to the oscillations of the bound state 
frequencies, the on-site superconducting order parameter 
at ri^2 also oscillates as the distance between the two im- 
purities is increased, as shown in Fig. EJJb). Note that a 
minimum in the oscillations of the order parameter co- 
incides with a larger splitting between the even and odd 
bound states. This is expected since a larger splitting 
implies a stronger coupling between the impurity bound 
states. At the same time, a larger coupling leads to a 
stronger effective scattering strength of the two-impurity 
system, and thus to a larger suppression of the supercon- 
ducting order parameter. 

With increasing scattering strength, JS^ of the two 
magnetic impurities, the bound states move to lower en- 
ergies, as follows from a comparison of the LDOS for 
JS/t = 2.0 [see FigCta)] with that for JS/t = 2.5 [see 
Fig.jS^a)], and eventually cross zero energy. For the case 
of a single magnetic impurity in an 5-wave superconduc- 
tor, it was first shown by Sakura:^^ that such a zero- 
energy crossing at a critical value, {JS)cr^ is the signa- 
ture of a first order quantum phase transition, in which 
the ground state of the superconductor changes from a 
state with spin polarization (Sz) = to a state with 
(Sz) = ±1/2, depending on the orientation of the mag- 
netic impurity. Here, the spin polarization is defined as 



LDOS [arb. units] 



SC order parameter [t] 



{Sz) = \jd\ p duj [N^{r,uj)-N^{r,co)]nF{oo) , 



where iVcr(r, cj) is the LDOS of the electrons with spin 
a =t, i, and ni? is the Fermi distribution function. This 
phase transition arises from a level crossing in the su- 
perconductor's free energy, JT, resulting in a discontinu- 
ity of dJ^/dJ at the transition; hence the first order na- 
ture of the transition (for a more detailed discussion see 
Ref.d^). For a single magnetic impurity, the transition 
also possesses a characteristic signature in the on-site or- 
der parameter which exhibits a discontinuous decreases 
and changes sign at {JS)cri thus undergoing a 7r-phase 
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FIG. 5: (Color online). LDOS at ri,2 as a function of energy 
(left column) and superconducting order parameter along r = 
(x, 0) (right column) for two impurities located at ri = (—2, 0) 
and rs = (2, 0) for (a)(b) JS/t = 2.5, (c),(d) JS/t = 2.7, and 
(e),(f) JS/t = 4.0. 



shift in comparison to its bulk value. 

For the case of two magnetic impurities, quantum in- 
terference leads to the emergence of two quantum phase 
transitions at {JS)cr,i and (J5')cr,2, corresponding to 
the transitions (Sz) =0 -1/2 and {Sz) = -1/2 ^ 
— 1, respectively, of the ground state's spin polarization. 
Specifically, the even bound state crosses zero energy 
at {JS)cr,i/t ^ 2.54, such that for JS > {JS)cr,i, the 
spin- 1 (spin- 1) peak of the even bound state is hole-like 
(particle-like), as shown in Fig.jSfc) for JS/t = 2.7. How- 
ever, while in the non-self-consistent T-matrix approach, 
the frequency of the even bound state, Cc;^ ^ evolves con- 
tinuously as a function of JS and reaches zero at the 
transition^^, within the BdG formalism, co^ ^ crosses zero- 
energy discontinuously at {JS)cr^ similar to the transi- 
tion associated with a single magnetic impurity^^ . Specif- 
(^4j ically, for {JS)cr,i — 0+ one finds ^^^^/t = 



while for {JS)cr,i + 0+ one has col^/t = ±0.001 [the 
upper (lower) sign corresponds to the spin-| (spin-|) 
peak]. Interestingly enough, we find that the energy of 
the odd bound state does not exhibit a discontinuity at 
{JS)cr,i' A similar discontinuous zero-energy crossing of 
the bound state energies is found at all phase transitions 
discussed below. The second quantum phase transition 
(with (Sz) = -1/2 -1) occurs at {JS)cr,2/t = 3.69 
when the odd bound state crosses zero energy, such that 
for JS > {JS)cr,2^ all spin- 1 (spin-|) peaks are now lo- 
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cated at positive (negative) energies [see Fig. |5Je) for 
JS = 4.M> {JS)cra]' 

In Figs. [3b), (d) and (f) we plot the superconduct- 
ing order parameter as a function of position for the 
three values of JS considered above. We find that the 
two phase transitions are accompanied by discontinuous 
changes of the on-site superconducting order parameter 
at ri^2, as follows directly from Fig. Efa) where we plot 
the on-site order parameter as a function of JS. Note, 
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FIG. 6: On-site superconducting order parameter as a func- 
tion of JS for two impurities located at (a) ri,2 = (t2, 0) with 
Ar = 4 (the critical values {JS)cr,j are indicated by arrows) 
and (b) n = (-1, 0) and r2 = (2, 0) with Ar = 3. 

however, that while the on-site order parameter is re- 
duced at {JS)cr,i^ it remains positive and does not ex- 
hibit a TT-phase shift, in qualitative contrast to the case of 
a single magnetic impurit y^^i^^i^^ . For JS > {JS)cr,i the 
superconducting order parameter decreases continuously 
with increasing JS, and crosses zero at JS/t ^ 3.15. 
This zero-crossing of the on-site order parameter does 
not lead to any signature in the LDOS. At (J5')cr,2, the 
superconducting order parameter at ri^2 exhibits a sec- 
ond discontinuity. Note, however, that in contrast to 
the discontinuous change at {JS)cr,i^ where the magni- 
tude of the superconducting order parameter decreases, 
its magnitude increases at {JS)cr,2- Interestingly enough, 
the suppression of the superconducting order parameter 
becomes spatially more confined when the system un- 
dergoes a phase transition. For example, the supercon- 
ducting order parameter at r^n = (3,0), i.e., one site 
away from one of the impurities, is much less reduced 
from its bulk value for JS > {JS)cr,i [Fig.[5fd)] than for 
JS < {JS)crA [Fig.ISi;b)]. Moreover, for JS > {JS)cr,2 
the order parameter at r^n as well as between the two 
impurities at r = (0, 0) is quite close to its bulk value [see 
Fig. [5ff)], in contrast to the results shown in Fig. [SJ^b) 

for JS < (J5)cr,l. 



For all interimpurity distances we considered, the de- 
pendence of the LDOS and the superconducting order 
parameter on the scattering strength is similar to the one 
discussed above. However, the splitting between {JS)cr,i 
and {JS)cr,2 depends strongly on the interimpurity dis- 
tance. For example, for Ar = 3, the two phase transi- 
tions occur at {JS)cr,i/t = 2.87 and {JS)cr,2/t = 2.9, as 
shown in Fig. |Sfb). The splitting between the two criti- 
cal values is thus considerably smaller than for the case 
Ar = 4 shown in Fig.|Hfa). We note, however, that while 
the splitting between the two critical values is small, they 
are nevertheless reduced from the critical value of a sin- 
gle magnetic impurity given by {JS)cr/t = 3.09. At the 
same time, we find that the energy splitting between the 
odd and even bound states is considerably smaller for 
Ar = 3 than for Ar = 4. In particular, for Ar = 3 
the frequency splitting between the even and odd bound 
state peaks at {JS)cr,i is Ao; ~ 0.0056t, while for Ar = 4 
at {JS)cr,i^ the frequency splitting is about one order of 
magnitude larger with Auj ^ 0.056t. We therefore con- 
clude that the splitting of the critical values of JS, similar 
to the splitting of the even and odd bound state energies 
depends on the effective coupling between the two im- 
purities. Hence, (J5')cr, 1,2 are functions of Ar. In this 
regard, it is interesting to note that with increasing JS, 
the frequency splitting between the odd and even bound 
states increases for Ar = 4, but decreases for Ar = 3, as 
can be seen from a comparison of Fig.^a) and Fig.[7fb). 
This effect directly refiects the dependence of the scat- 
tering phase shift on JS, as discussed above. 

In Fig. d^a) we present the superconducting order pa- 
rameter at ri^2 for several values of JS as a function of 
inter-impurity distance, Ar. For JS/t = 2, the order pa- 
rameter exhibits only weak Friedel-like oscillations when 
the distance between the impurities is varied (see also 
Fig. EJ. In contrast, for JS/t = 2.7, the on-site super- 
conducting order parameter, A^ 2? oscillates much more 
strongly when the interimpurity distance is changed, for 
example, from A^ 2 = 0.030t at Ar = 3 to A^ 2 = 0.003t 
at Ar = 4. The reason for this strong variation becomes 
clear when one considers Fig. For Ar = 3, one finds 
{JS)cr,i/t > JS/t = 2.7, implying that the spin polariza- 
tion of the superconductor is (Sz) =0. In contrast, for 
Ar = 4 one has {JS)cr,i/t < JS/t = 2.7 < {JS)cr,2/t, 
yielding (S^) = -1/2. In other words, for JS/t = 2.7 
the superconductor undergoes a first order phase tran- 
sition when the interimpurity distance is changed from 
Ar = 3 with spin polarization (Sz) = to Ar = 4 where 
(Sz) = — 1/2. As a result, the superconducting order 
parameter varies strongly between Ar = 3 and Ar = 4, 
as shown in Fig. [Tfa). This effect directly refiects the 



dependence of {JS)c 



•,1(2) 



on Ar, as discussed above. 



While for JS/t = 2.7, the spin polarization of the su- 
perconductor's ground state changes between {Sz) = 
and (Sz) = —1/2 , we find that for JS/t = 3.0, its spin 
polarization varies between (Sz) = — 1/2 [for Ar = 2, 4, 6 
where {JS)cr,i/t < JS/t = 3.0 < {JS)cr,2/t] and (Sz) = 
-1 [for Ar = 1,3,5 where {JS)cr,2/t < JS/t = 3.0]. 
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FIG. 7: (Color online), (a) Superconducting order parameter, 
An 2 ^1,2 as a function of interimpurity distance, Ar, for 
several values of JS/t. (b) Frequency of the even and odd 
bound state peaks as a function of inter- impurity distance 
Ar for JS/t = 2.5. 



Finally, for JS/t = 4.0, the behavior of the supercon- 
ducting order parameter is complementary to that at 
JS/t = 2.5, in that the superconductor's spin polariza- 
tion is (5'^) = -1/2 for Ar = 2, and (S^) = -1 for 
Ar 7^ 2 (while (S^) = for JS/t = 2.5 and Ar ^ 2). 

In order to explore the interplay between the position 
of the bound state peaks in the LDOS and the supercon- 
ducting order parameter further, we plot in Fig.[7fb) the 
energies of the even and odd bound states as a function 
of inter-impurity distance for JS/t = 2.5. For Ar = 2, 
the spin polarization of the superconductor is given by 
{Sz) = —1/2, since the odd bound state has crossed zero 
energy. At the same time, the superconducting order 
parameter at ri^2 is significantly reduced [see Fig. Q^a)]. 
In contrast, for all other interimpurity distances, no zero- 
energy crossing of the bound state peaks is observed, and 
the superconducting order parameter at ri^2 is signifi- 
cantly larger than that for Ar = 2. Both results imply 
that for Ar ^ 2, one has (5'^) = 0. 



IV. THREE MAGNETIC IMPURITIES 

We next consider the case of three magnetic impuri- 
ties, for which the superconducting order parameter ex- 
hibits an interesting behavior at the first order phase 
transitions. Specifically, we study three impurities with 
parallel spins (S||z) that are located at ri^a = {^x^O) 
and T2 = (0,0). Quantum interference of scattered elec- 
trons again lifts the degeneracy of the fermionic bound 
states associated with each of the impurities, leading to 
the emergence of six peaks in the LDOS. In order to 
gain insight into the symmetry of the hybridized bound 
states, we assume that the unhybridized bound states of 
energy Eq are coupled to the nearest neighbor (next near- 
est neighbor) state via a hopping term K {K'). The re- 
sulting energies of the three hybridized bound states are 
then given by 

Eodd = Eq — K' 

The wavefunction with energy Eodd is given by ^odd = 
{ipi — ips)/y/2 where ipi is the wavefunction of the unhy- 
bridized fermionic bound state induced by the impurity 
at r^. Note that the wavefunction "^odd is odd and thus 
vanishes at any point in space with equal distance to 
ri and ra. Thus the odd bound state vanishes in par- 
ticular at r2, which, as discussed below, is important 
for the behavior of the superconducting order parame- 
ter at the first order phase transitions. The wavefunc- 
tions with energies E± for the case = are given by 
= \/2{tlji ± i\)2/\p2. — ?/^3)/>/5, respectively (we omit 
their forms for K' since these are quite cumbersome 
and irrelevant for the discussion below). Note that the 
energy of the odd bound state always lies between those 
of the ^±-states, with the energy difference between the 
^±-states and the ^o^^^^-state given by 

A^± = ^± - EoM = ^ ± y ^ + 2/^2 . (5) 

Next, we consider the concrete case of three impurities 
located at ri^a = (t2,0) and r2 = (0,0) with nearest 
neighbor distance Ar = 2. The resulting LDOS as a 
function of frequency at ri^3 for four different values of 
JS is shown in Fig. [SI As discussed above, one finds 
six peaks in the LDOS, three peaks each arising from the 
components of the hybridized bound state wave functions 
with Sz = +1/2 (spin-t peaks) and Sz = —1/2 (spin-| 
peaks). For small JS [Fig.ISfa)], all spin-| (spin-|) peaks 
are hole-like (particle-like) and thus located at > 
{co < 0). With increasing JS, the peaks move towards 
zero energy, and a phase transition occurs when the first 
bound state crosses zero energy at {JS)cr,i/t = 1.98. At 
this first order transition, the spin polarization of the 
superconductors changes from (Sz) = to (Sz) = —1/2. 
For JS = 2M > {JS)cr,i [see Fig.E^b)] one of the spin-i 
(spin- 1) peaks is located at negative (positive) energies. 
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FIG. 8: (Color online). LDOS as a function of frequency at 
ri,3 for the case of three impurities with interimpurity dis- 
tance Ar = 2 and four different values of JS\ (a) JSfi — 1.9, 
(b) JSjt = 2.5, (c) JSjt = 4.0, and (d) JSjt = 4.5. The 
results in (a) were obtained for a lattice with — 34 and 
Ny — 22, those in (b)-(d) for a lattice with Nx — 28 and 
Nx = 18 (see Ref.^^). 



The spatial pattern of the three bound states at the 
frequencies of the spin-| and spin-| peaks is shown in 
Fig. El The bound state whose spatial pattern is shown 
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FIG. 9: (Color online). Spatial intensity plot of the LDOS 
for three impurities with JSjt — 2.5 at the frequency of the 
spin-| peaks and spin-| peaks of the ^--state [(a) and (b)], 
the ^odd-state [(c) and (d)], and the ^+-state [(e) and (f)], 
respectively. 

in Fig. I^c) and (d) is easily identified as the odd bound 
state with wavefunction ^ odd- In order to identify the 
^±-bound states, we use Eq.(|S|) and the results presented 



in Figs. IHa) and[7|;b). Since K = D{Ar = 2) > 0, 
and K' = D{Ar = 4) < 0, we obtain from Eq.(|S|) 
that |A£;_| > |A£;+|. An inspection of the LDOS in 
Fig. IHl then allows us to identify the bound state whose 
spatial pattern is shown in Figs. [HJa) and (b) as the 
^_-state [which crosses zero energy at (J5')cr,i], while 
the bound state shown in Figs, ^e) and (f) is the 
state. This result is also supported by our analysis of 
the LDOS-- based on the value of ^± at ri,2. Fi- 

nally, we note that the spatial structure of these bound 
states remains almost unchanged as JS is increased and 
the superconductor undergoes a first order transition. 

The zero energy crossing of the ^_ -state at {JS)cr,i is 
accompanied by a discontinuous jump in the supercon- 
ducting order parameter at all three impurity sites, as 
shown in Fig. UnT a). While the superconducting order 
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FIG. 10: (Color online). On-site superconducting order pa- 
rameter at ri,2,3 as a function of JS for three impurities 
with nearest neighbor distance (a) Ar = 2 (the critical values 
{JS)crj are indicated by arrows) and (b) Ar = 3. 

parameter is reduced at the transition, it remains posi- 
tive, and does not undergo a 7r-phase shift. When JS is 
further increased, the odd bound state crosses zero en- 
ergy at {JS)cr,2/t = 3.72, and the superconductor un- 
dergoes a first order phase transition to a state with 
(5'^) = — 1 [a spatial intensity plot of the LDOS at fre- 
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quency of the odd bound state is shown in Figs-IH^c) and 
(d)]. For JS = 4.0t > (J5')cr,2, two spin-| (spin-|) peak 
are located at < (cc; > 0) [see Fig. IHI^c)]. Since the 
odd bound state vanishes at r2, we find that the super- 
conducting order parameter at r2 remains unaffected by 
the zero energy crossing, while that at ri^s changes dis- 
continuously and undergoes a 7r-phase shift, as shown in 
Fig.CUa). Finally, at {JS)cr,3/t = 4.31 the ^+-bound 
state crosses zero energy, and the superconductor under- 
goes a first order phase transition to a state with a spin 
polarization (Sz) = —3/2 [the spatial LDOS intensity of 
the ^+ bound state is shown in Figs. Elje) and (f)]. All 
spin-| (spin-|) peaks are now located at a; < (a; > 0), as 
shown in Fig.I^Jc) for JS = AM > {JS)cr,3- At (J5')cr,3, 
the superconducting order parameter at all three impu- 
rity sites again changes discontinuously, and increases in 
magnitude, as shown in Fig. ITnT a). However, only the 
order parameter at r2 changes sign at this transition. 

The behavior of the superconducting order parameter 
at the first order phase transitions is different from the 
above scenario when the nearest neighbor interimpurity 
distance is increased to Ar = 3 and the impurities are 
located at ri^a = (^3, 0) and T2 = (0, 0). In this case, the 
three phase transitions occur in a much smaller range of 
JS than for Ar = 2, with {JS)cr,i/t = 2.73, {JS)cr,2/t = 
2.79 and {JS)cr,3/t = 3.01. As discussed above, this 
indicates a weaker coupling between the unhybridized 
bound states of the impurities, and we therefore also ex- 
pect a smaller frequency splitting between the peaks of 
the three bound states in the LDOS. This smaller fre- 
quency splitting is directly evident from a comparison of 
the LDOS for Ar = 3 and JS/t = 2.78 > {JS)cr,i shown 
in Fig. [ni with the LDOS for Ar = 2 shown in Fig. |H1 
Moreover, the fact that the energy splittings between the 
peaks, |A£^±|, is almost identical, implies that the effect 
of K' = D{Ar = 6) can be neglected in comparison to 
that of = D{Ar = 3). A more detailed analysis of 
^± then shows that for Ar = 3, the relative energies of 
the ^+ and ^_ bound states are reversed in comparison 
to the case Ar = 2: the peaks of the ^_ bound state 
are now located at higher energies than those of the ^+ 
bound state. Thus the ^+ bound state crosses zero en- 
ergy at (JS')cr,i, the ^odd bound state at (J6')cr,2, and 
the ^_ bound state at {JS)cr,3- Interestingly enough, 
the LDOS shown in Fig. [TTT a) at ri^a only possesses five 
peaks. The reason for the "missing" peak becomes ev- 
ident when one considers the spatial form of the LDOS 
at the bound state frequencies. Specifically, the LDOS of 
the ^+ bound state at the frequency of the spin-t peak 
vanishes at ri^a, as shown in Fig. [TTT c). which explains 
the absence of the sixth bound state peak in Fig. [TTT a). 
Note that the LDOS at the spin-| peak of the ^+ bound 
state does not vanish at vi^s [see Fig. lTTT b)]. This spatial 
form of the ^+ bound state possesses an interesting con- 
sequence for the behavior of the superconducting order 
parameter at (J5')cr,i, as shown in Fig. ITnT b). Specif- 
ically, at {JS)cr,i the superconducting order parameter 
at T2 is reduced, while the order parameter at vi^s is 
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FIG. 11: (Color online). LDOS for three impurities with 
Ar = 3 and JS/t = 2.78 > {JS)cr,i- (a) LDOS at ri,3 as 
a function of frequency. Spatial intensity plot of the LDOS 
at the frequency of the spin-| peaks and spin-j peaks of the 
^+-state [(b) and (c)], the ^odd-state [(d) and (e)], and the 
^--state [(f) and (g)], respectively. 



slightly enhanced. This behavior of the superconducting 
order parameter is quite different from that at (J5')cr,3, 
when the ^_ bound state crosses zero energy. Since nei- 
ther the LDOS of its spin-| peak nor of its spin-| peak 
vanishes at vi^s (or r2), we find that at {JS)cr,3^ the order 
parameter at all three impurity sites decreases discontin- 
uously. Moreover, when the LDOS of both the spin-| 
and spin- 1 peak vanishes at an impurity position, as is 
the case with "^odd bound state at r2 [see Figs. lTTT d) and 
(e)], the superconducting order parameter at this impu- 
rity site is unaffected by the zero energy crossing of the 
bound state, as follows from Figs, a) and (b). Thus 
the vanishing of the LDOS of only one of the spin-t or 
spin- 1 peaks at an impurity site can lead to a discontin- 
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uous increase in the superconducting order parameter at 
the phase transition, as discussed above. 

V. CONCLUSIONS 

In conclusion, we have studied quantum interference 
phenomena in impurity structures consisting of two or 
three magnetic impurities located on the surface of an 
5-wave superconductor. By using a self-consistent BdG 
formalism.^^, we simultaneously investigated how quan- 
tum interference affects (a) the formation of impurity 
bound states in the LDOS, and (b) the spatial form of 
the superconducting order parameter. 

We obtained three important results. First, we find 
that the characteristic signatures of quantum interference 
in the LDOS and in the superconducting order param- 
eter are coupled. In particular, the larger the splitting 
between the hybridized bound states in the LDOS, the 
stronger is the suppression of the superconducting order 
parameter at the impurity sites. As a result, both the 
frequency of the bound states as well as the on-site or- 
der parameter exhibit Friedel-like oscillations as the dis- 
tance between impurities is changed. Second, the sup- 
pression of the superconducting order parameter around 
impurities does not qualitatively change the signatures of 
quantum interference in the LDOS, thus demonstrating 
the robustness of quantum interference phenomena. The 
physical origin of this robustness lies in the fact that, 
while the superconducting order parameter is suppressed 
by impurities, it recovers its bulk value over a length 
scale which is set by the Fermi wavelength, A^, and not 
by the superconducting coherence length, which is typi- 
cally much larger than A^. Third, we show that quantum 
interference leads to multiple first order quantum phase 
transitions in the superconductor, which are accompa- 
nied by a change in the spin polarization of the supercon- 
ductor's ground state. The superconductor can be tuned 
through these transitions by increasing the impurities' 



scattering strength or by changing the interimpurity dis- 
tance, the latter being the experimentally more relevant 
tuning. These quantum transitions exhibit several char- 
acteristic features that are qualitatively different from the 
ones of a phase transition associated with a single mag- 
netic impurity. In particular, the superconducting order 
parameter does not necessarily decrease discontinuously 
or undergo a 7r-phase shift, and, in certain cases, can even 
be enhanced at the transition. The difference between 
the critical values of JS for which these multiple transi- 
tions occur is determined by the hybridization strength, 
and is thus directly linked to the splitting between the 
impurity bound state energies. The tuning of the super- 
conductor's spin polarization by changing the interimpu- 
rity distance potentially possesses important applications 
in the field of spin electronics^^ and quantum information 
technology^^ as it might lead to the creation of new types 
of quantum qubits. 

Finally, we note that the relative orientation of the im- 
purities' magnetic moments is determined in general by 
the residual interactions between the impurities, whose 
consideration is beyond the scope of this study. In case 
where the dominant contribution to the inter-impurity 
interaction arises from an RKKY-type process, the in- 
teresting question arises to what extent the interaction 
is affected by the formation of the hybridized impurity 
states. This question will be addressed in some future 
work. 
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